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Abstract. We present a method for calculating explicit expressions of the shear three-point function for various 
cosmological models. Fhe method is applied here to the one-halo model in case of power law density profiles for 
which results are detailed. The three-point functions are found to reproduce to a large extent patterns in the shear 
correlations obtained in numerical simulations and may serve as a guideline to implement optimized methods for 
detecting the shear three-point function. In principle, the general method presented here can also be applied for 
other models of matter correlation. 



1. Introduction 
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f^*) . Cosmic shear surveys offer a unique and unbiased view on the way matter cluster at cosmological scales and can then 
serve as a mean to explore the details of the gravitational instabilities. In particular the statistical properties of the 
O ■ shear field, in particular their high order correlation functions, are expected to reflect those of the matter density field 
(see Mellier 1999). The measurement of correlation functions in the cosmic shear surveys has then been recognized as 
a mean to break the parameter degeneracy between the density parameter of the Universe and og that measurements 
^ ' of the shear power spectra provide for (Bernardeau et al. 1997, Jain & Seljak 1997) or a way to test the gravity law 
d ■ at cosmological scales (Bernardeau 2004). 

Matter correlation functions have been explored extensively in different regime. A reasonable understanding of their 
• *h . behaviors have been obtained in the quasilinear regime from perturbation theory techniques and into the nonlinear 
K% regime from phenomenological approaches (see review of Bernardeau et al. 2002). Recently a semianalytic technique 
'■ based on the halo model (Seljak 2000, Ma & Fry 2000, Scoccimarro et al. 2001, see Cooray & Sheth 2002 for a 
recent review) has attracted a lot of attention because it can provide accurate predictions for the galaxy or matter 
correlation properties. So far however most of these theoretical approaches have dealt with the matter density field 
or with the convergence field. The later is closely related to the density field, being a simple projection of it, but the 
actual observational situation is more cumbersome. What is directly measured is the shear field (through galaxy shape 
measurements, see Mellier 1999 for a detailed review on the observational techniques) and it reveals difficult to build 
reliable projected density maps. It forces us to explore more deeply the high-order shear correlation patterns. Those, 
because of their geometrical nature, are more intricate than that of the convergence maps. In particular complex 
relations are expected between the correlation components because both of pseudo- vector nature of those and because 
one expects vanishing "B" modes in the shear field. As a results shear correlation coefficients depend in a complicated 
way on the basis on which they are computed, have frequent sign changes, etc. (see Schneider and Lombardi 2003, for 
an exhaustive exploration of these properties) . It makes the construction of methods for measuring the shear three- 
point function difficult. Attempts however have been made to detect the shear or aperture mass high order correlation 
functions (Bernardeau et al. 2003, Pen et al. 2003, Jarvis et al. 2004). Although these methods can give robust results, 
they may be far from optimal and therefore not completely satisfactory. One aim of this paper is to provide better 
foundations for the design of new methods for measuring shear three-point functions. 

In this paper the focus will be put on the one-halo model case for the description of the matter correlation 
properties. The idea is not however to build an accurate model for the mass distribution but to build a model good 
enough to reproduce the generic geometrical properties of the shear correlation patterns. The 1-halo model, with a 
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power law density profile, is a good candidate for such a task being simple enough to allow complete computations 
of the correlation functions. The paper is divided as follows. The second section is devoted to a presentation of the 
mathematical tools applied to the one-halo model in the third section. The fourth section gives some insights on the 
results that have been obtained. 



2. A complex variable representation 

What we are interested in are the statistical properties of the either the projected density field <f>(r), the convergence 
field n(r) and the 2-component shear field (71,72). In the following the plane wave approximation is used, as usual, so 
that the convergence and shear field are all obtained through derivative operators applied to the projected potential, 

«(r) = {dl + dl)^{v) (1) 
(71,72) - {dl-8l2d a d b ) <£(r) (2) 

where r is the angular position of the point on the sky and a and b are the point coordinates. 

In principle the convergence or shear correlations are functions of the point coordinates. Given the statistical 
isotropy and homogeneity of the fields those correlation functions should depend only on the relative positions and 
are globally rotationally invariant. Those correlation functions can therefore be written for instance as a function of 
the three distances of equivalently as a function of 2 distances and one angle, etc. The coordinate system advocated 
in this paper, and the one I will exploit throughout, is however none of those. Let me denote the point coordinates as 
ri, r 2 and r 3 respectively ; l\, I2 and l 3 are the three lengths of the triangle they form such that Z§ = |r 2 — i"i| 2 , etc. 
Then I can define the two ratios, 

;2 /2 



72' » J ;2 
f 3 *3 



(3) 

and the complex variable 

z 3 = ^ (l - x 3 + y 3 - -x 3 + y 3 ) 2 - 4y 3 ) . (4) 

It should be noted that (1 — x 3 + y 3 ) 2 — 4y 3 < so that z 3 has generically a non vanishing imaginary part. The idea 
is then to write all correlation functions as a function of (k,Zi), i being either 1, 2 or 3. It is first to be remarked that 
Zi is nothing but the components of r, in the basis defined by the other two points 

Z3 -» («3 -ib 3 )/h- (5) 

As a result the derivatives with respect to point coordinates can simply be re-expressed as derivatives with respect to 
Zi and Zi: 

dri = hd Zt +d St ,-id iei +id St ) (6) 

H 

It implies that the Laplacian operator, or operators that give access to the shear coordinates out of the potential field, 
respectively read 

dl + d 2 b = ^d z ck (7) 
dl-dl = (8) 

2d a d b = -p(^ 2 -a|). (9) 

Clearly those equations suggest that the shear field is better written as a complex number 7 = 71 + i72 and we then 
have 

1 

P 

-^d z 2<j>(li,Zi) and 7^,2-) = - 

Full computations with such variables will however been tractable only if change coordinates can be written in a 
natural way. Indeed the transform rules that give the Z\ and z 2 , complex coordinates of r x and r 2 , as a function of z 3 
are actually simple to compute and read, 

21 = 7^— . z 2 = -^— = l--. (12) 
1 - z 3 1 - zi z 3 



K{k,Zi) = pdzidzMk, Zi) (10) 
l{k, Zi) = -Ld z 2(p(li,Zi) and 7(^1, z») = -4 % *») ( n ) 
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The coordinate change rules have to be completed with the following relations, 

^1=^3 ZiZi, l2=h z 2Z2 (13) 

that gives the transformation rules for the base lengths. These rules suffice for the computations of correlation of scalar 
quantities, e.g. correlations between the projected potential and the convergence fields. 

For the computation of correlations between shear components however, further transform rules are required. These 
are the ones that rotate the shear coordinates from one basis to another. After a bit of vector analysis, these changes 
of coordinates simply amount to multiply the complex shear coordinates by (1 — z 3 )/(l — z 3 ) in order to change base 
(ri,r2) to (r2,r 3 ) or by similar operations for the other changes (applying circular permutations to the indexes). 

If the projected potential correlation function is known, then the computation of either the convergence field or 
the shear field is straightforward. What it takes is only a series of proper derivatives and rotation operators applied 
to the potential correlation function. 

This approach is particularly attractive when the latter correlation function can be written as a sum of terms of 
the form f(z)f(z). In this case indeed the operators act independently on each factor. As we will see next this is 
precisely the situation that is encountered in the one-halo model with power law profiles. 



3. Correlation functions in the 1-Halo model 

The halo model provides with a simplified description of the matter correlation properties. It has been advocated to 
be an accurate approach for describing the galaxy and matter correlation properties and studies of the dark matter 
density field through halo contributions have recently been developed in Seljak (2000), Ma & Fry (2000), Scoccimarro 
et al. (2001), but also applied to three-point lensing statistics (Cooray et al. 2000b and Takada & Jain 2003). 

In this section we will illustrate the methods sketched in the previous section for shear three-point functions using 
the halo model restricting for simplicity to the one-halo term. The latter is indeed expected to dominate on the 
small scales (see for instance Fosalba et al. 2005) where the three-point function of the shear is easiest to measure in 
observations. In itself the use of the halo model for the computation of the lensing high order correlations is not new. 
It has been done by Cooray & Hu (2001) for the convergence bispectrum and in Takada & Jain 2002a, Zaldarriaga 
& Scoccimarro (2003) for the shear correlation functions. The accuracy of the model and the has been extensively 
commented in those papers and in Fosalba et al . (2005). It won't be commented any further it here. 

The basic relation for computing the correlation functions in the one-halo model has been provided for by 
Scoccimarro (1997, taking advantage of a result given by Davydychev 1992) as discussed in Zaldarriaga and Scoccimarro 
(2003). 

The one-halo amounts to assume that a stochastic field <j> can be represented by a randomly located halo with a 
spherically symmetric </> profile, <p(r) ~ F(\r — ro|) where ro, the center of the halo is evenly distributed in space. The 
4> correlation functions then read, 

F„(n...r„) = (#-i)... #r n ))~ J d d ro^(|r 1 -r |)...F(|r„-r |) (14) 

This expression is in general complicated. We are interested in such an expression for the three-point function 
and in case of a power law profile, F(x) ~ x~ a . It is to be noted that such an integral converge for d — 3a < and 
d — a > 0. For F(x) = F x~ a the final result can be given in a closed form in terms of hypergeometric functions. This 
result has been presented by Scoccimarro (1997) for the three-point function taken advantage of integrals developed 
in the context of high energy physics. 

We are interested in this result in the context of a 2D projected density field 1 . The result reads 

7T 3 

F 3 (n,r 2 ,r 3 ) = l\ 3a - sin(mr/2)2 cos(mr/2) r(2 _ 3a/2 )r( a /2)2 X (15) 

[I>/2)r(l - 3a/2)G(z 3 ) G(z 3 ) - T(2 - 3a/2)r(l - a/2)H(z 3 )H(z 3 )} (16) 
where the variables are those defined in the previous part and where G{z) and H{z) are respectively, 

G{z) _(!-*) — (17) 

TJ( ~\ „l-a 2Fi(f,l-f,2-Q!,z) ( 

g(z) = z r(2-a) (18) 

1 The scales at which the correlations are to be measured are small enough so that the small angle approximation can be 
used. 
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This expression results from a slight modification of the formula presented in Scoccimarro (1997). From the form 
<|16|) it is clear that F 3 is real. 

The idea pursued in this paper is that convergence fields can be described to a reasonable accuracy with the 1-halo 
model with a density profile close to a = 1. The goal though is not to use (fltSf) for the convergence field but to construct 
in this framework the shear three-point functions, 



£*'ft7k( r i> r 2> r 3) = (7i(ri)7i( r i)7fe( r i)) 



(19) 



where i,j, k stand for the two component indices of each of the shear field. 

To obtain this function two routes can be followed. The one advocated in Zaldarriaga and Scoccimarro (2003) was 
to transform the formal expression of £ 7i7j . 7fc into a sum of integral each of which can be computed in a manner similar 
to F 3 . It leads however to expressions with a very large number of terms that reveal hard to manipulate and simplify. 

The route followed here is radically different. As suggested in the second section it amounts to take proper sub- 
sequent derivatives of the potential three-point function, ^^(ri, r2, r 3 ) with respect to the positions. The potential 
correlation function is a priori not well defined because it requires the use of (| 1 fc>l) when a is about -1. It demands then 
a regularization of the profile at large distance. It can be checked however that the introduction of this cutoff does not 
affect the finite distance behavior of its derivatives if those are finite. 

The general mathematical operators to be used have been presented in the previous section. Least by not last, 
the calculations remain tractable in case of this model thanks to the fact that the function G and H have simple 
derivatives. It can indeed be shown that (using Gradshteyn and Ryzhik 1980, Eqs. 9.137.6, 9137.7 and 9.137.17) 



d n 
dz n 



G(z) 



d n 

dz n K ' 



G {n \z) = {l-z) 
H {n \z) 



i- a -n 2F1 (f , 1 - f , a + n, z) r(a/2 + n) T(l - a/2 + n) 

T{a/2) T(l - a/2) 



1 9- 
1— a— n * 



1 



r(a + n) 
|r, 2 — a — n, z) 



r(2-a-n) 

To get the shear correlation functions the subsequent steps should be followed, 
1 



i<l><fo{h, z 3 ) 

z 2 ) 

£777(^1' z \) 



,2 d z^£,<t><t><t>( l 3,Z 3 ) 
l 3 

J2^z 2 {R^t^ih, z 2 )} 

'2 

-^13 72^1 {-^21^77(^)^1)} 
60 



(20) 
(21) 

(22) 
(23) 
(24) 



and similar expression for the other components of the shear correlation functions, £7(1-1)702)7^3) i £7(ri)7(r 2 )7(r 3 )) etc.. 
Note that in those expressions a change of variable should be applied from one line to another (for both the z and 
the I variables) ; the rotation operators are R21 = I/R12 = (1 — z\)/(l — z i), R13 = I/-R31 = (1 — z 3 )/(l — z 3 ), 
-R32 = -R31-R12 = (1 — z 3 )/(l — z 3 ) (1 — zi)/(l — zi) and they should be replaced by their inverse (or equivalently their 
complex conjugate) each time it applies to conjugate components of the shear. 

After calculations and simplification taking advantage of recursion relations of the hypergeometric functions and 
with the help of formal calculators one can obtain explicit expressions for the shear correlation functions. The first 
four read 



?7(ri)7(r 2 )7(r 



C a {h)a 2 {2 + a) 2 



{-l + z 3 fzl 



{ 2 a (1 - 2 z 3 ) 2 (-2 - z 3 + zi) + 8 (l - 2z 3 + 



Az% 



2zl) + 



3 a' -4- 



3 zi- 16 zt 



*4)] r(|)r(-i- r ^)G( )(z 3 )G (0) ^3)+ 



- a) 
3 a 



■ 5z 3 

2 - 5 z 3 + 2 zj + 2 4 - 5 zl + 2 z\) V(±) Y{-1 + ^) G« (z 3 ) G (0) 



8(1 

r(2-^)r(l-|) [(2a(l-2z 3 ) 2 (-2- 
2 (-4 + 5z 3 + 3z 2 



3 a N 



23 



(1-2 



Z3 



-4z 3 3 -2z 3 4 ) 



3 a 
8(1 



-4- 
a) (2 



164 + 8z 3 4 )) H^(z 3 )- 



5z 3 



2z\ 



2z\ 



5z* 



2zl) HU(z 3 ) 



H (o) 



^7(ri)7(r 2 )7(r 3 ) 



C a (l 3 )a 2 (2 + a) 



r(~) r(-i + ^) ((-2 + 3 a) (8 (-1 + z 3 ) z 3 + 



(-1 + z 3 ) z 3 (-1 + 23) z 3 
a (-I-8Z3 + 8Z3 2 )) G (0) (z 3 ) + 4 (-1 + 2 23) (2 (-I + Z3) z 3 + 



(25) 



Vy(ri)7(r 2 )7(r 3 ) — 
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a (-l-2z 3 + 2z 3 2 )) G«(z 3 )) ((-2 + 3 a) G (0) (l 3 ) + 4 (-1 + 2 z 3 ) G«(z 3 ))" 
r(2-^)r(l-|) [(-2 + 3a) (8 (-l+z 3 ) z 3 + a (-1 - 8 z 3 + 8 z 3 2 )) ff<°>(* 3 )+ 
4 (-l + 2z 3 ) (2 (-l + z 3 ) z 3 + a (-1 - 2 z 3 + 2 z 2 )) ff (1) (z 3 )] ((-2 + 3 a) # (0) (z 3 )+ 
4 (-l + 2« 3 ) # (1) (^))} 

C a (Z,)a(2 + a) / 2 r( l + |) F(-l + ^) [(-2 + 3a) (-4 z 3 + a (-1 - 5 z 3 + 2 z 3 2 )) 
(-I + 23J z 3 z 3 L z z 

G (0) (« 3 ) + 4(1 + z 3 ) (-2 z 3 + a (l - 4 z 3 + z 3 2 )) G (1) (z 3 ) 

(-2 + 3a) (l + 2z 3 ) G<°)(l 3 )+4 (-l + z 3 2 ) G«(z 3 ) 
ar(2-^)T(l-|) [(-2 + 3 a) (-4z 3 + a (-1 - 5 z 3 + 2 z 2 )) H(°\z 3 )+ 
4(1 + z 3 ) (-2 z 3 + a (1 - 4 z 3 + z 2 )) ffW (z 3 ) x 

(-2 + 3a) (l + 2z 3 ) i/(°)(z 3 ) + 4 (-l + * 3 2 ) H^(z 3 )]] 

C " (/3)Q(2 + Q) 3 ( 2 r d + ?) r M + ^) f(-6 + « (9 - 6^3) + 4z 3 ) G(«)(z 3 )- 



3a, , a , 

ar(2-— )r(i--)x 



with 



C a (/ 3 ) 



(-l + z 3 ) (-I+T3) z 3 6 I 2' 2 

4 (-2 + z 3 ) z 3 G (1) (z 3 )] [(-2 + 3 a) (4 (-l + z 3 )+a (-4 + z 3 + 2 z 3 2 )) G<°>(Z 3 ) + 
4 (-2 + z 3 ) (2 (-l + z 3 ) + a (-2 + 2z 3 + z 3 2 )) G (1) (z 3 ) 

(-2 + 3a) (-3 + 2z 3 ) ( ) (z 3 ) + 4 (-2 + z 3 ) z 3 i?< 1 )(z 3 ) 

(-2 + 3 a) (4 (-l + z 3 )+a (-4 + z 3 + 2 z 3 2 )) H^(z 3 )+ 

: (-2 + z 3 ) (2 (-l + z 3 )+a (-2 + 2z 3 + z 3 2 )) ff (1) (z 3 )]} 



(tt 3 csc(^) 2 sec(^) 



2T(2- ^)r(f)' 



(26) 



(27) 



(28) 



(29) 



The other four components can be derived as complex conjugates. The results are given as a function of z 3 and z 3 and 
therefore when the coordinates (in particular the shear coordinates) are expressed in the (ri,r2) basis. 



4. Shear Patterns in the One-Halo Model 

The amplitude of the shear correlation function can then be compared to the amplitude of the convergence correlation 
function. To start with one can consider £7(2) defined as 



1/2 



ri)tj (r 2 )7fc(r 3 ) S7i (ri)7j (r 2 )7fc (r 3 ) 

ijk 



(30) 



It is to be noted that there are some symmetry properties in this quantity due to the fact it should be invariant 
under z — ► 1/(1 — z) and z — > 1 — 1/z when properly rescaled (these changes correspond to transform that leaves the 
triangle shape unchanged; the amplitude of the shear functions are then only changed by a scale factor effect). 

Its amplitude in shown on Fig. ^ where the spatial dependence of the amplitude of the shear correlation function 
is compared to that of the convergence correlation function. Whereas the convergence correlation function has a very 
simple spatial dependence, the spatial dependence of the shear correlation is much more complex exhibiting minima 
for certain configurations. It is to be noted however that these patterns are quite dependent of the power spectrum 
index (right panel of Fig. 1 and Fig. 2). In particular the existence of a local maximum of the signal for equilateral 
triangles vanishes for steep spectra. 
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-0.6-0.3 0. 0.3 0.6 0.9 1.2 1.5 -0.6-0.3 0. 0.3 0.6 0.9 1.2 1.5 



Fig. 1. Comparison of the spatial dependence of the amplitude of the three point function for both the convergence 
fied (left panel) and the shear field (right panel). These result correspond to the one-halo model with a = —1.. In this 
plot, ri = (0,0) and r 2 = (1,0). As can be observed the spatial structure of the shear three-point function is much 
more complex than the one of the convergence field 




Fig. 2. Comparison of the spatial dependence of the amplitude of the three point function for a = —0.75 (left panel) 
and a = —0.5 (right panel). 



A fair fraction of the signal however can be extracted from a subpart of the correlation functions. Following 
Bernardeau et al. (2003) it is possible to map at least some aspects of the shear patterns while considering the 
function, 

6(ra) = {(7^-7^)7^) (31) 
which can be written in terms of the complex coordinate correlations as 

6( r 3) = 2 (^7(ri)7(r2)7(r3) + %(r 1 )7(r 2 ) 7 (r 3 )) ( 32 ) 

The £3 patterns are depicted on Fig. EI They are very similar to what has been observed in numerical simulations. 
This part of the components does not however encompass the signal coming from equilateral configurations. The 
extension of the uniform and strong pattern in between the two points depends again on the halo profile. 
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-0.6 -0.3 0. 0.3 0.6 0.9 1.2 1.5 -0.6 -0.3 0. 0.3 0.6 0.9 1.2 1.5 



Fig. 3. Pattern of the £3 field in case of the one-halo model for Z3 — 1, a — —1. (left panel) and a = —0.5 (right panel). 
In this plot, ri = (0, 0) and r2 = (1,0). They level indicates the amplitude of this component of the shear correlation. 

5. Prospects 

The use of complex coordinates proved fruitful in computations of shear correlation functions. It has been possible to 
obtain the shear three-point correlation functions in manageable forms for the one-halo model with a power law profile. 
In the context of this model those results can be used to infer general properties for the shear correlation patterns, as 
in Benabed & Scoccimarro (2005); a first step toward the construction of robust methods for the measurement of the 
shear three-point function. 

The theoretical framework presented in this paper is actually very general and can be used for alternative de- 
scription of the density correlation properties as such those used in the quasilinear or the nonlinear regime. Indeed 
the mathematical relation between the projected potential correlation functions and those of the projected density is 
rather simple and could be integrated out for simple cases 2 opening the way for explicit computations of the shear 
correlation functions in more diverse cases. 
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